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Abstract. 

Based on the Liouville-Weyl definition of the fractional derivative, a new direct 
fractional generalization of higher order derivatives is presented. It is shown, that 
the Riesz and Feller derivatives are special cases of this approach. 
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1. Introduction 

The fractional calculus PQ, [4] provides a set of axioms and methods to extend the 
coordinate and corresponding derivative definitions in a reasonable way from integer 
order n to arbitrary order a: 

Qn Qa 

The definition of the fractional order derivative is not unique, several definitions e.g. 
the Riemann, Caputo, Liouville, Weyl, Riesz, Feller, Grunwald fractional derivative 
definition coexist [7|-[T2]. In the last decade, there has been a steadily increasing 
interest in applications of the fractional calculus on such different fields of research 
like mechanics, anomalous diffusion or fractional wave equations. 

Most work is dedicated to the special case of the first order derivative operator 
(n — 1) replaced by an appropriately chosen fractional derivative operator. This 
approach is tempting in that sense, that higher order derivatives my be replaced in 
a natural way by a consecutive sequence of first order derivatives and consequently 
may be replaced by the corresponding sequence of the fractional extension of the first 
order derivative. 

Until now, there exists no general approach for a direct fractional extension of 
higher order derivatives, except for (n — 2). For that case, Riesz and Feller have 
derived a fractional generalization of the second order derivative operator directly. 

In this letter, we will derive a general definition of a direct fractional generalization 
of higher order derivatives. For that purpose, we will first collect the necessary tools, 
which are currently used for a fractional generalization of the first- and second order 
derivative. We will then propose an extension of the fractional derivative definition to 
higher order derivatives. 
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The Liouville-Weyl fractional integrals of order < a < 1 are defined as 

iu(x) = _L /" (a _ o a -vm ( 2 ) 



r(a) 

r(a) A 

The Liouville-Weyl fractional derivatives of order < a < 1 are defined as the left- 
inverse operators of the corresponding Liouville-Weyl fractional integrals 



D%<j>(x) = I+ a cj>(x) = +Aji-^ (a;) (4) 

D«<t>{x) = lZ a m = -JLll-°<f,(x) (5) 
The definitions (jl} and (J5]) may be written in an alternative form: 

' T(l - a) J c q+ 

" n ; r(l-a) Jo r +1 
which may be derived via 

D%cj>{x) = I^ a <j>(x) = + ( 8 ) 



1 d 
F(l -a)dx J_ 

i a 



T(l - a) y 



(x-o _a 0(Ode (9) 



r Q <M*-£R (io) 



1 r r(4^-0K (ii) 



T(l-a)J * " dC 



r(i - a) u e +l y e 

A specific linear combination of the Liouville-Weyl fractional integrals results in the 
Riesz fractional integral J^: 

Ir<Kx) = + 7 7^ 0*) = / l*-ei a_1 ^(0de a >0, a ^1,3(51) 
2cos(a7r/2) 7_ oc 

The Riesz fractional derivative is then given by 

D a R et>{x) = - ^t D L ^) (14) 
2cos(a7r/2) 

or, according to (JB]),®: 

sin(c»r/2) [°° #s + fl - 20(g) + 0(s - fl n ^ 
fflW^) = r(l+a) y i^pj d£, < a < 2(15) 

Feller has proposed a generalization of the Riesz fractional derivative of the form 

I$<t>(x) = (c_(0,a)I? +c+(0,a)I« )<t>{x) (16) 
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sm(™) 

c+(0 ' a) = — iSP) — (18) 

The Feller fractional derivative is defined as 

D^s) - - (c+(0, a)Dl + c-{9, a)D a _) 4>{x) (19) 
Setting 9 = we obtain 

c_(0 = O,a) = c + (0 = O,a) = - ^— (20) 

2cos(a7r/2) 

which coincides with the definition of the Riesz fractional derivative (fl4|) . 
Another special case results for setting 6 = 1 

C -(fl = l,«) = -c + (fl=l,a)= , 1 (21) 

2 sm(o!7r/2) 

which leads to the simple form of the fractional derivative: 

D?Hx) = D +~°^ Mx) (22) 
2 sin(a7r/2) 



= T(l + a) ' I -dH (23) 



( cos(a7r/2) f°° <f>(x + - <p{x - f ) 
T° 

This derivative should be interpreted as the regularized Liouville-Weyl fractional 
derivative (fT4| . 

Therefore the Feller fractional derivative may be rewritten as a linear combination 
of Df and D%: 

D$4>(x) = (Ai(0, a)(D% - D a _) + A 2 (9, a)(D« + D°)) <j>(x) (24) 

with 

At (8, a) = - \ (c+ (9, a) - c_ (8, a)) = - - , / sin(07r/2) (25) 

1 2sm(a7r/2) 

A 2 (9,a)= - \ (c + (9,a)+c.(9,a)) =- ^~J- /9 , cos(fl7r/2) (26) 

2 2cos(air/2) 

which finally reads: 

D%<j>(x) = (sm^Tr/ 2 )- ? + cos(9n/2)D%) <f>{x) (27) 

In this form of the Feller fractional derivative the parameter 9 may be interpreted 
rather as a rotation parameter instead of a skewness parameter, proposed by other 
authors. In addition, this form is better suited for a generalization to higher order 
derivatives, which will be performed in the following section. 

3. The fractional generalization of higher order derivatives 

Using the basic properties of the central differences operators 

6 i tKx)=<t>(x + ±t)-<Kx-±t) (28) 

h<Kx) =^(cj>(x + 0-<f>(x-0) (29) 
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we define the central differences operator J) fe of order k 



for k even 
for k odd 



or explicitely, using (|28p and (f^^jl : 

2[(fe+l)/2] 

D fc 0(x)= ^ a^-([(fc + l)/2]-n)0 

n=0 

with the summation coefficients 

k 
n 

1 |Y fc-1 \ _ / fc-1 

2 I n ) ~ { n-2 



a « = (-1) 



for k even 
for k odd 



The renormalized fractional derivative is then given as: 

1 



D*' a <t>(x) 

and the normalization factor 
'p+i)/2] 



ir»V(*) 



aw e + 



T(l-a) 



J2 a k n (k-n-l) a 



n=0 



cos(7ra/2) 
sin(7ra/2) 



for k even 
for k odd 



(30) 



(31) 



(32) 



(33) 



(34) 



With (|33[) based on the Liouville definition of the fractional derivative we therefore 
have given all fractional derivatives, which extend the ordinary derivative of order k: 

d k 

lim k D a = — 

a^k dX K 

In addition, for these derivatives the invariance of the scalar product follows: 



( k D a *f*(x))g(x)dx = (±) h 



f(x)*( k D a g(x))dx 



The first four fractional derivative definitions according (f3"3"|) follow as: 



in>)=r(1 + Q) ^) x 



f( x + Z)-f(x-Z) 



(35) 
(36) 

(37) 



a+l 



„, \ „, , sin(a7r/2) 
2 D a f(x) = r(l + a) i '—L x 



< a < 1 



(38) 



/(a; + £)-2/(z) +/(*-£) 



a+l 



3 ^/(,) = r(i + a) cos(OT/2) . 1 , 



< a < 2 



7T 2 Q - 2 
/(a: + 20 + 2/(z + - 2f(x - $ + f(x - £) 



a+l 



(39) 



< a < 3 
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„ „, „ , % sin(a7r/2) 1 
4 D a f(x)=T(l + a) ^ 7 ; 2a _ 4 x (40) 

3 -/(a + 2fl + 4/(s + Q- 6/(g) + 4/(x - fl - /(x - 2Q 

< a < 4 

(41) 

These definitions are valid for < a < &. Setting a > k 

jn/c 

fci?a = 5^*^""* < neN (42) 

and chosing n so that < a — nk < k the definitions given are valid for all a > 0. 

In the same manner the Feller fractional derivative definition may extended to 
fractional derivatives of higher order. 

We introduce hyper sperical coordinates on the unit sphere on R n : 

xi = cos(0„_i) (43) 
x 2 = sin(6>„_i)cos(6>„_ 2 ) (44) 

%-i= sin(0„_ 1 )sin(0 r( _2)...cos(0i) (45) 
x n = sin(0 n _i)sin(0 n _ 2 )...sin(0i) (46) 

With these coordinates the Feller definition of a fractional derivative may be extended 
to 

n 

*D<* {ek} = Y,XkkD a (47) 
fc=i 

4. Conclusion 

Based on central differences a generalized fractional derivative of arbitrary order has 
been propsed. 
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